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The one- dimensional propagation of waves in a bichromatic potential may be modeled by the 
Aubry- Andre Hamiltonian. The latter presents a delocalization-localization transition, which has 
been observed in recent experiments using ultracold atoms or light. It is shown here that, in contrast 
to Anderson localization, this transition has a classical origin, namely the localization mechanism is 
not due to a quantum suppression of a classically allowed transport process. Explicit comparisons 
with the Anderson model, as well as with experiments, are done. 
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I. INTRODUCTION 

It is by now well established that disorder may lead 
to the absence of diffusion in the propagation of linear 
waves (Anderson localization). The phenomenology of 
Anderson localization strongly depends on dimensional- 
ity. In one dimension (ID) and for uncorrelated disorder, 
it has been shown that all wave functions are exponen- 
tially localized. This implies that the transmission of 
a wavepacket incident on a disordered region of length 
L decreases as ex.-p{—L/Lioc)^ where Lioc <C L is the 
localization length. This phenomenon is valid for arbi- 
trary energies, even when the energy E of the incident 
wavepacket is much higher than the typical amplitude 
of the disordered potential. In this respect, Anderson 
localization is a strongly non-classical phenomenon. 

Another class of interesting potentials are quasiperi- 
odic potentials. They arise, typically, from the superpo- 
sition of two periodic lattices with incommensurate peri- 
ods. Part of their interest is that they provide an exam- 
ple of intermediate potentials between ordered (e.g. pe- 
riodic) and disordered potentials: the system has neither 
translational symmetry nor true disorder. In particular, 
the model we consider here (the Aubry-Andre Hamil- 
tonian) is known to present a delocalization-localization 
transition as a parameter is varied; wavefunctions change 
from extended to exponentially localized 

In recent experiments both random [2] and quasiperi- 
odic potentials [sj have been used to directly observe 
the localization of matter waves. The expansion of 
non- interacting Bose-Einstein condensates through one- 
dimensional quasiperiodic potentials was shown to clearly 
display a delocalized- localized crossover (with a shift, 
however, with respect to the predicted critical value of 
the parameter) Q. A similar transition has been re- 
cently observed with light propagating in quasiperiodic 
photonic lattices [i]. Our purpose here is to provide a 
semiclassical analysis of the experiments and, more gen- 
erally, of the Aubry-Andre Hamiltonian, which allows 
for a simple and lucid picture of the underlying physi- 
cal phenomena (for related semiclassical analysis see e.g. 
[H, IBI, 0, H S 0])- In particular, and in contrast to 
what is quite often claimed, the analysis shows that there 



are deep fundamental differences between the Anderson 
localization and the localization in quasiperiodic poten- 
tials, although both present exponential localization of 
the wavefunctions. 

As already pointed out above, Anderson localization is 
a strongly non-classical phenomenon, where a quantum 
particle of energy E incident on a ID disordered region 
of length L ^ oo is reflected with probability one, even 
though its incident energy is much higher than the maxi- 
mum height of the disordered potential. Thus, in Ander- 
son localization a quantum particle is perfectly reflected, 
whereas its classical counterpart is perfectly transmitted. 
In this sense, Anderson localization is a purely wave- 
mechanical effect. In contrast, as shown below, localiza- 
tion in quasiperiodic potentials is different. It is a purely 
classical effect, not due to any interference or wave-like 
dynamics. Indeed, already at the classical level the tra- 
jectories display a delocalized-localized crossover. In this 
integrable model the wavefunctions can directly be asso- 
ciated to some specific classical orbits, around which they 
are exponentially concentrated. As a parameter is var- 
ied, they simply follow the classically driven transition. 
Quantum mechanically, there is thus no suppression of 
a classical diffusion process. It is our purpose to stress 
these aspects of the model, and to clarify the connections 
with recent experiments. 

II. THE MODEL 

The Hamiltonian we are interested in is of the tight- 
binding form, with nearest-neighbor hopping between 
discrete position sites Qn and a superimposed quasiperi- 
odic potential 

H J'^{\qn+i){qn\^\qn){qn+i\) 

n Q\ 

+ WY,co^{2^f3n)\qn){qn\. 

n 

where jS is the commensur ability parameter, \qn) is a lo- 
calized Wannier state on site n of the lattice, and J and 
W are parameters that control the amplitude of the hop- 
ping transitions and of the quasiperiodic potential, re- 
spectively. Since the Hamiltonian is invariant under the 
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increase of f3 by an arbitrary integer, we can restrict to 
< /3 < 1. This Hamiltonian is known as the Aubry- 
Andre or Harper Hamiltonian in the litterature [l|, [l2[ • 
There are two well known physical problems effectively 
described by Eq. ([T]). The first one is the motion of elec- 
trons in two dimensions for a periodic potential in the 
presence of a magnetic field applied perpendicular to the 
plane, when interband transitions are neglected [13]. In 
this case P is related to the ratio of the area of the mag- 
netic flux quantum to the unit cell in coordinate space. 
The second one, directly related to the experiments de- 
scribed here, is the ID motion of particles in the presence 
of two superimposed periodic potentials, the main one of 
period Ai that determines the position of the discrete lat- 
tice points ^n, and the perturbing one of period A2; in 
this case /3 = A1/A2. The connection between the latter 
problem and Eq. ([T]) is explicitly described in the next 
section and in Appendix A. 

Expressing the hopping between nearest-neighbor sites 
of the main lattice in terms of the translation operator 
ex.-p{—iXip/h) {q and p are the usual position and momen- 
tum operators, respectively), one can transform Eq. ([T]) 
into the more symmetric form 



H = 2Jcos(27rj9/P) + ly cos(27r(7/Q) , 



(2) 



where P = 27rfi/Ai and Q = X2. The classical analog of 
this Hamiltonian is obtained by replacing the operators 
by c-numbers 

n{q, p) = 2 J cos(27rp/P) + W cos{27Tq/Q) . (3) 

This classical Hamiltonian is periodic in both position 
and momentum. Therefore its study can be restricted to 
a single cell of size (Q,P). Note that 7i{q^p) contains, 
through the momentum scale factor P = /i/Ai, quantum 
Planck's constant. Its validity and applicability, as well 
as the apparent paradox that it contains a quantum scale, 
will be discussed below. 

TL{q^p) defines a time-independent one-dimensional 
problem which conserves the energy. It is thus an in- 
tegrable system, with Eq. ([3|) defining, for different val- 
ues of the energy E = H{q^p)^ one dimensional curves 
in the phase space plane (q^p). The classical phase- 
space dynamics is sketched in Fig. [1] for different val- 
ues of a = W/J. For a = all trajectories are ex- 
tended in space (e.g., the projection of the classical curves 
onto the q axis covers the entire interval [0,(5], in con- 
trast to localized orbits that cover only a fraction of 
it), and have constant momentum, whereas for a 1 
exactly the opposite happens, trajectories are localized 
in space and delocalized in momentum. In between, 
a crossover is observed, with the presence of two sep- 
aratrices that delimit three different types of trajecto- 
ries. The first type are closed localized orbits that os- 
cillate around (q^p) = ((5/2, P/2). This is a standard 
clockwise oscillation around the minimum of cos{27rq/Q). 
The second one are also closed localized oscillations, now 
around {q^p) = (0,0), mod((5,P), which is a less stan- 
dard counterclockwise oscillation around the maximum 
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FIG. 1: Phase space trajectories of Hamiltonian (|3]) for differ- 
ent values of the parameter a = W/ J (full lines). In the Tight- 
est column, the Husimi representation of several eigenstates 
for /3 = 1/51 is superimposed to the classical trajectories (see 
text). 



of cos(27rg'/(5) due to a local negative mass. Finally, the 
third type, in between the separatrices, are extended or- 
bits, in g for a < 2 and in p for a > 2. At a = 2 the two 
separatrices merge together. This is the critical value of 
the parameter a above which all classical trajectories are 
localized in space. 

The quantum mechanical motion is richer than the 
classical one, because of the additional parameter /3 = 
A1/A2 = 27rh/QP. Its presence in the quantum Hamil- 
tonian is a consequence of the discrete nature of the mo- 
tion in the q direction, defined by the position of the 
main lattice sites, and the commensur ability effects with 
respect to the secondary lattice. The quantum oper- 
ator Eq. ([2]) commutes with the translation operators 
Tq = exp{—iQp/h) and Tp = exp{iPq/h)^ which trans- 
late by one unit cell in each phase-space direction. How- 
ever, a quantum dynamics simultaneously periodic in 
both directions, like the classical one, is possible if and 
only if the translation operators also commute with each 
other. The most general condition enforces commutation 
by translations over a phase-space domain of dimensions 
(M(5,P), where M is an arbitrary positive integer. It is 
easy to show [9] that the operators Tmq and Tp commute 
if and only if is a rational number 
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(4) 



where N is an arbitrary positive integer. Having in mind 
the old Weyl rule "one quantum state per phase-space 
volume /i", the physical interpretation of this condition 
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is clear: in order to be periodic in both directions, the 
extended phase-space cell (MQ^P), made of M units in 
the q direction, should contain an integer number N of 
quantum states. This defines the torus {MQ^P) as the 
quantum phase space. An arbitrary quantum state now 
satisfies the generalized boundary conditions 



TMQm = e 



(5) 



where = (6>i,6>2) are good quantum numbers pre- 
served by the dynamics (the Bloch phases). The Hilbert 
space hence breaks down into N dimensional subspaces 
parametrized by ^. H is therefore a finite-dimensional 
periodic operator whose spectrum is absolutely continu- 
ous. The latter consists of N Bloch bands EiiO)^ defined 

by 



Q,...,N-l 



(6) 



where are the eigenvectors that satisfy Eqs. ([5]). 

Strictly speaking, due to the periodicity, all eigenstates 
are thus delocalized for any value of the parameter a. 

Note that the semiclassical limit is obtained when 
N ^ oo^ i.e. when the number of quantum states is large 
(with, possibly, M ^ oo but M/N ^0). In this case the 
fraction of phase space occupied by one quantum state 
tends to zero, and the quantum dynamics is described 
with increasing precision by the classical Hamiltonian, 
Eq. (|3]). The semiclassical limit is, thus, equivalent to 
the limit ^ (or, more generally, and due to the pe- 
riodicity, to the limit /3 tending to an arbitrary integer 
from above). In practice, /3 plays the role of an effec- 
tive Planck's constant. Since /3 = A1/A2, an appropriate 
choice of the period of the two superimposed periodic 
potentials allows to tune /3 and control how quantum 
or classical the system is. This limit is thus different 
from the usual semiclassical limit, where one formally 
takes h ^ 0, or considers a high energy regime. From 
this point of view, there is thus no contradiction in the 
fact that the classical Hamiltonian H{q,p) contains /i, 
since the path towards the classical behavior described 
by that Hamiltonian is controlled by the two frequen- 
cies Ai and A2. This emphasizes that the classical limit 
P ^ we are considering has a more general character. 
It is one example among a wider class of semiclassical 
limits, sometimes generically called large- A/" limits, that 
are often encountered in physics (see for instance [14] and 
references therein). In our case, the classical limit does 
not correspond to the standard limit of a classical parti- 
cle of kinetic energy /2m moving in the presence of two 
superimposed periodic potential (cf Eq. (p!Q|) below). Al- 
though purely classical, the structure of H{q,p) keeps in- 
formation about the original quantum mechanical model, 
Eq. ([1]). The periodicity in the p direction arises from the 
discretness of the lattice sites distant by Ai, imposed by 
the tight-binding form of the dynamics related to a tun- 
neling process. The boundedness of the kinetic energy 
term has important consequences, in particular in the 



appearance of the unusual localized in space-delocalized 
in momentum classical trajectories observed in the limit 
of a large a. 

Now we turn to the case of irrational whose quan- 
tum mechanical behavior deserves special attention. In 
contrast to our previous discussion, in this case the phase 
space is not compact, the Bloch angles are not good quan- 
tum numbers, and the matrix H is of infinite dimension. 
For any irrational value of p it can be shown [l5| that 
the spectrum is a Cantor set (among which the famous 
Hofstadter butterfiy for a = 2). In addition, the model 
displays a delocalization- localization transition [l|, 
For Of < 2, the spectrum is absolutely continuous and all 
states are extended. For a > 2, the spectrum is pure 
point and all states are localized. Finally, at a = 2 the 
spectrum is singular continuous, with multifractal eigen- 
states [1]. 

The ideal situation concerning any physical experiment 
whose aim is to display the delocalization-localization 
transition is to implement an irrational /3, typically the 
golden mean (>/5 — l)/2. In practice, one can only ap- 
proach an irrational as some rational approximation. It 
is thus of interest to analyze an irrational P as the limit 
of a sequence of rational numbers. The most efficient 
sequence (in terms of convergence) is known to be the 
continuous fraction expansion 



mi 



(7) 
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This sequence of rational approximations Pi = 1/mi, 
P2 = 777,2/ (mi m2 + 1), • • • generates a sequence of pe- 
riodic systems that approximate the quasiperiodic one 
with increasing accuracy. In the lowest order approxima- 
tion Pi = 1/mi, only one unit cell is quantized (M = 1), 
with N = mi quantum states supported by this cell. For 
such a situation, semiclassically the WKB method allows, 
in ID, to construct the eigenstates tlJi{0) and eigenvalues 
Ei{0) from the quantization of some specific classical tra- 
jectories, in a one-to-one correspondence. The mi states 
are selected by the Bohr-Sommerfeld quantization rule 
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S= I pdq = (3QP [ i + 



^2 



(8) 



where i = 0, . . . , mi, (wi^ w^) are the winding numbers 
of the quantized classical trajectory t in the (g,p) di- 
rections, respectively, and S — S{E) is the action of t. 
The Maslov index 1/2 applies only to trajectories with 
{wi^W2) = (0,0) (closed trajectories). The eigenvalues 
of states associated to (0, 0) trajectories are thus, at this 
semiclassical level of approximation, independent of 0. 

Eq. ([8]) associates one quantum state to one classical 
trajectory of the elementary phase space cell, and de- 
fines the spectrum of energies Ei{0). In this descrip- 
tion, the associated quantum state '0i(6>) can be shown 
to strongly concentrate or localize around the corre- 
sponding quantized classical trajectory. To illustrate 



this point, it is useful to display the eigenstates in the 
phase-space Husimi representation [9], defined in Ap- 
pendix B. In this representation, to each eigenstate V^i(^) 
of energy Ei{0) corresponds a positive-definite phase- 
space function Wi{q^p; 9) which, in the semiclassical limit 
= mi oo, behaves as pjj 
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state 200 



Wi{q,p;e)^ 



v{q,p) 



exp 



j_ / n{q,p)-Ei{e) 

fni I v{q,p) 



(9) 

where v{q^p) = ^ fP" is the phase-space velocity of 
the classical trajectory of energy Ei. The rightest column 
of Fig. [U superimposes to the classical trajectories the 
corresponding quantized eigenstates, for = mi =51, 
M = 1, ^ = (0,0) and different values of a. One can 
clearly see, for a given a, the concentration over the cor- 
responding classical trajectory as well as, for varying a, 
the transition from extended to localized states, which 
simply follow the classical crossover. For open trajecto- 
ries (either in the q direction for a < 2 or in p for a > 2) 
the quantum state is supported by two symmetric and 
isoenergetic classical curves, coupled by tunneling. This 
effect is observed at ^ = (0, 0), where the eigenstates are 
real. The quantum state combines both classical trajec- 
tories, each of which has a current, but whose overall 
superimposed net current vanishes. As d is varied, the 
eigenstate become complex and can concentrate on one 
or the other of the classical trajectories. 

At the second step of the rational approximation to 
an irrational /3 ^ (32 = 7x12/ {mim2 + 1), the quan- 
tum phase-space is an extended torus of dimensions 
{m2Q,P). This space now supports N = mim2 + 1 
quantum states. With respect to the previous approx- 
imation /3 ^ /3i = 1/mi, the size of the quantum cell 
and the number of states N ^ mim2 increase by a factor 
m2] in contrast the number of quantum states per unit 
phase-space cell remains almost constant. The scheme 
repeats as one increases further the order of the con- 
tinuous fraction approximation. Along this process, the 
size of the extended phase-space cell increases in the q 
direction (thus tending towards a cylindrical geometry), 
the dimension of the corresponding Hilbert space also 
increases, while the effective Planck constant which 
measures the fraction of the elementary phase space cell 
occupied by one quantum state, remains approximately 
constant. At a given step of the approximation and a 
given ^, from Eq. (|8|) it can be shown that two types of 
quantum states exist, classified according to the struc- 
ture of the wavefunction in the extended phase space 
(MQ^P). First, the delocalized ones, which are asso- 
ciated to classical orbits open in the q direction of the 
(±1,0) type (that exist only for a < 2). Second, lo- 
calized states, associated to either closed orbits (of the 
(0, 0) type) or to open orbits in the p direction (of the 
(0,±1) type). These states are typically localized over 
one (or several, when resonances occur) isoenergetic clas- 
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FIG. 2: Husimi functions for N=431 and M=516 for two ar- 
bitrary quantum states as a function of a. Their evolution 
clearly illustrates the delocalized- localized crossover. 



sical orbits of the extended phase space. Their width, 
controlled by remains almost constant along the con- 
tinuous fraction approximation. Thus, as the size of the 
extended torus increases in the q direction, the relative 
width of the localized states, compared to MQ, dimin- 
ishes. Husimi plots of wavefunctions, that confirm our 
general picture, can be found in Ref. [To( for sequences 
of rational approximations of the golden mean. 



III. RELATION TO EXPERIMENTS 

In recent experiments [3|] , a cloud of non-interacting ul- 
tracold ^^K atoms created from a Bose-Einstein conden- 
sate is released in a bichromatic potential produced by 
superimposing two lasers of wavelength Ai and A2 (and 
of corresponding wavenumbers ki and /C2) in a standing 
wave configuration. The dynamics is described by the 
following Hamiltonian 

i^B = - — ^ + Vi cos^(fcix) + V2 cos^{k2x) , (10) 

where Vi and V2 denote the amplitude of each of the 
two laser beams and m is the mass of the atom. If the 
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primary lattice (say, Vi) is deep enough, that is if the 
recoil energy Eji = h'^kf/2m is much lower than Vi, one 
can approximate the motion in the primary lattice by a 
tunneling process between localized states of neighboring 
sites, each of which feels the presence of the superimposed 
additional periodic potential V2. The Hamiltonian can 
thus be reduced to the tight-binding form Eq. ([1]), where 
/3 = A1/A2. The effective parameters J and W in Eq. ([1]) 
are directly related to Vi, V2 and Er (see Appendix A). 

Several different laser frequencies have been used in 
the experiments. In [3] they have realized the system 
described by the Hamiltonian (p!Q|) with Ai = 1032 nm 
and A2 = 862 nm, which yields A1/A2 = 516/431 (431 
is a prime number). Restricting to < < 1, we have 
/3 = 85/431 0.1972. Moreover, the typical amplitude 
of the primary lattice is ten times larger than the corre- 
sponding recoil energy, which makes ([T]) a good approxi- 
mation to the experiment. With the above values of Ai, 
A2 and the quantum system under study is therefore 
periodic in both q and p directions (because /3 is rational), 
with = 431 states accomodated in M = 85 elementary 
cells. This means that every 431 wells of the primary 
lattice, the system repeats. The Hamiltonian ([1]) is thus 
a finite matrix of dimension A^ = 431, and the eigen- 
values and eigenstates depend parametrically on the two 
phases ^ = (^1,^2)- It can be shown (see Ref. [9] and 
Appendix B) that the two phases are determine by the 
relative position of the two lattice sites. Though not di- 
rectly measured in the previous experiments, the Bloch 
phases are fixed by the experimental conditions. To be 
specific, we fix them here to ^ = (0,0). Coming back to 

its continuous fraction decomposition gives two lower 
order approximants. Pi = 1/5, and P2 = 14/71. The sys- 
tem is thus nearly periodic after five wells of the primary 
lattice. The approximation p ^ Pi means also that ap- 
proximately five quantum states are accomodated in one 
elementary phase space cell {Q^P). Moreover, the rel- 
atively low value of p ~ 0.2 indicates that the classical 
approximation to the dynamics, given by Eq. (|3|), as well 
as the the semiclassical picture developped in the pre- 
vious section, is a meaningful framework, as discussed 
below. 

Fig. [2] shows the Husimi function of several eigen- 
states of Eq. ^ represented in the quantum phase space 
(MQ^P) for the experimental conditions, namely p = 
85/431, and different values of a, assuming 6 = (0,0). 
One verifies that the eigenstates follow, for different a's, 
the general description given at the end of last section. 
For a <C 2 the eigenstates are concentrated along the 
delocalized classical trajectories (mixing right and left 
propagating orbits, because 6 = (0,0)), while for a 2 
the eigenstates are concentrated on localized trajecto- 
ries that wind around in the p direction. In between a 
crossover between these two limits is observed. This fig- 
ure illustrates, under the experimental conditions, that 
the transition from delocalized to localized states as a 
increases may be interpreted in terms of classical dynam- 
ics, and is therefore different in nature from the Anderson 
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FIG. 3: RMS width of a cloud of ^^K atoms after t=750 ms 
of expansion in a (quasi) periodic potential characterized by 
p = 85/431 = 0.1972 ... as a function of a. The initial width 
of the cloud is 5 /xm, indicated by a horizontal dashed line. 
Dots : experimental results for different values of J/h (taken 
from Ref. [3]). Dashed-dot : theory (see text). 



localization. Moreover, as already pointed out in the in- 
troduction, this behaviour is generic and is not restricted 
to the special set of parameters used in the experiment. 

Experimentally, the transition was probed through the 
evolution, for different values of a, of an initial cloud 
of atoms released from the ground state of a confin- 
ing magnetic trap. The initial state is thus a Gaussian 
wave packet centered around some initial point qo and 
with zero average initial momentum, with an extension 
of about ten main lattice sites. After letting the cloud 
expand for some time (750 ms), the final RMS width of 
the wave packet is measured. Their results are reported 
in figure [31 To test our setting, we have computed, using 
the eigenstates discussed above, the time evolution |^(t)) 
of such an initial state, given by 

N-l 

= exp ^-iE,{0)t/h) ^,{zo)\MO)) , (11) 

where, for the actual experiment, A^ = 431, V^^(zo) is the 
coherent state representation of the eigenstate i, where 
^0 = Qo/V^ (cf Appendix B), and the Bloch angles 6 
are set to zero. In Fig. [3] the width averaged over ini- 
tial position computed from Eq. (pT]) is compared to the 
experimental data. The agreement is quite good. Note 
that the transition takes place here for a ^ 2, as pre- 
dicted by theory, and not around a 7, as obtained in 
Ref. [3]. The difference is due to an improved estimate 
of the dependence of the effective parameters W and J 
upon Vi, V2 and Er (see appendix A and ^]). Note also 
that the maximal experimental width observed after ex- 
pansion, of the order of 40/im, corresponds to ~80 main 
lattice sites, much smaller than the phase space cell, that 
covers 431 main lattice sites. 
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IV. COMPARISON WITH THE ANDERSON 
MODEL 

In the previous sections we have stressed the classical 
origin of the delocalization-localization transition, and of 
the localized phase, in bichromatic potentials. Let us 
be more explicit, in particular in the comparison with re- 
spect to the Anderson model of localization, by emphasiz- 
ing their differences and similarities (differences between 
the quantum motion in disordered and quasiperiodic po- 
tentials in terms of effective dispersion relations have 
been recently analyzed in Ref. [l9|]; see also Ref . llQll V 

The original Anderson model of localization [201] with 
diagonal disorder is described, in ID, by a Hamitonian 
similar to Eq. ([T]) . The difference is that the energy of site 
n, instead of being cos(27r/5n), is now given by a random 
variable Vn having, for instance, a uniform distribution 
between -1 and 1. Then W simply controls the ampli- 
tude of the random potential, and WVn has a uniform 
distribution between —W and W. 

If one whishes to follow a semiclassical analysis simi- 
lar to the previous one but for the Anderson model, the 
first difficulty is to define appropriately the classical limit. 
The potential, as defined above, is not continuous, and 
has no direct classical limit. One way to get round this 
difficulty is to smooth the previous random potential, by 
for instance convoluting Vn with a Gaussian of width a 



V{q) = W^V„ 



e-(9-9n)/2<T" 



V2' 



(12) 



where Qn denotes, as in previous sections, the position of 
the sites of the discrete lattice separated by a distance 
Ai. Then, the classical limit of the Anderson model is 
described by the Hamiltonian 



n{q,p) = 2Jcos(27rp/P) + WV{q) , 



(13) 



where P = h/Xi. The classical phase space trajectories of 
such a Hamiltonian are, qualitatively, not very different 
from those of the quasiperiodic potential Eq. (|3]), in par- 
ticular in the two extreme limits a ^ and a ^ oo (see 
Fig. [4]): they are given, again, by delocalized constant-p 
trajectories in the former limit, and localized, constant-g, 
trajectories in the latter one (we are using, as before, the 
parameter a = W/ J) . The main difference is in the inter- 
mediate regime, where trajectories are more irregularly 
shaped and a larger number of separatrices is observed. 
Thus, as in the quasiperiodic Hamiltonian, in the Ander- 
son model a delocalized to localized crossover is observed 
classically as a increases. 

As we have seen in previous sections, the quasiperiodic 
case qualitatively follows the classical behavior: a quan- 
tum crossover from a delocalized regime to a localized 
one is observed as a increases. This contrasts with the 
Anderson model, where states are localized for arbitrary 
values of a > 0. Fig. 2] illustrates this point, by showing 
several quantum states for different values of a. We have 
numerically diagonalized the Anderson model on a finite 
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FIG. 4: Husimi distributions (gray code) of several eigenstates 
and their corresponding classical orbits (solid lines) in phase 
space for the Anderson model with N — 401 and a — \i. As 
a is increased the classical model also exhibits a delocalized- 
localized crossover. Localized quantum eigenstates already 
exist for a < 2, even if the corresponding classical orbits are 
delocalized in space, as shown for a = 1. 



chain of = 401 sites with = (0, 0) and obtained the 
classical phase space orbits from (p!3|) width a = Xi. For 
the sake of clarity Fig. [4] does not show all the classical 
orbits at a given energy but only the ones directly as- 
sociated to a quantum state. In the two limiting cases, 
namely a = and a ^ 1 (here 100), the quantum states 
are associated to very simple classical orbits, i.e. fully 
delocalized and fully localized trajectories, respectively. 
However, at finite a the eigenstates are localized, even 
those whose energy correspond to delocalized classical 
orbits, as shown in Fig. [H for a = 1 (of course, the length 
of the system should be longer than the corresponding 
localization length). 

The main difference between the quantum behavior of 
the two models is thus in the range a < 2, where the An- 
derson model, in contrast to the classical behavior, has 
localized states, while the quasiperiodic model follows the 
classical trajectories. For large values of a, the local- 
ization in the Anderson model has, as the quasiperiodic 
potential, a classical interpretation, because the classical 
orbits are themselves localized. In this respect, the non- 
trivial regime in the Anderson model is the non-classical 
localization observed for small values of a, an effect which 
is absent in the bichromatic potential. This stresses the 
differences and similarities between the two models, in 
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their quantum and classical behavior. 



and thus 



V. CONCLUSION 

In the present paper we have analyzed, using semi- 
classical methods, the delocalized-localized transition ob- 
served in quasiperiodic potentials. A general analysis of 
the corresponding Aubry- Andre Hamiltonian in terms of 
the commensurability properties of the parameter /3 was 
given. Particular attention was devoted to the interpre- 
tation of recent cold atom experiments, as well as to a 
comparison with the one-dimensional Anderson model of 
localization. 

Several extensions and generalizations of the dynamics 
considered here are of interest, with possible experimen- 
tal realizations. They are motivated by the possibility to 
explore more complex transport effects. Among them, we 
can mention the possibility to pulsate in time one of the 
two superimposed potentials, thus introducing chaotic 
motion in the dynamics described by the kicked Harper 
Hamiltonian [21]. Along similar lines of research, it is 
worthwhile to mention recent cold atom experiments con- 
cerning the observation of the Anderson metal-insulator 
transition [22[. Another important effect is related to 
the role and modifications induced by interactions in the 
delocalization-localization properties in quasiperiodic, as 
well as random, potentials [3 [13 • 

This work was supported by grants ANR-05-Nano- 
008-02 and ANR-NT05-2-42103 and by the IFRAF In- 
stitute. 



APPENDIX A: CORRECTION TO THE 
HAMILTONIAN PARAMETERS 

In Ref. \S\ they use the following notations. The orig- 
inal Hamiltonian describing the motion of cold atoms in 
a bichromatic potential is written 

Hb = ~~^~Q^ + SiEr^ COS^(/Cix) + S2Er^ COS^(/C2x) , 

(Al) 

where ki = 27r/A^ and Er. = fi^/c|/2m, z = 1,2, denote 
the wave numbers and the two recoil energies, respec- 
tively. Its reduction to a tight-binding form is expressed 
as 

H = J ^ \Wm^i){Wm \ + \Wm){Wm+l\ 
m 

+ W^cos(2T:l3m)\wm){wm\, 



(A2) 



where (3 = X1/X2 (not restricted to [0,1]). The parame- 
ters in (|A2p . namely J and W, were calculated as follows 

J = 1.43s?-^^exp(-2.07Vir) , 



W 



(A3) 



a = 



W 



S2 



2.86 



exp (2.07\/sr) 



(A4) 



On the other hand, in Ref. [2?] the mapping between 
the two Hamiltonians (|Aip and (|A2p have been studied 
very carefully, in order to analyze the presence of possible 
mobility edges due to deviations from the tight binding 
approximation. Writing their formulae with the same set 
of notations than in Ref. [3] yields 



W_ 
T 



8 .V^ 



exp (2 



(A5) 



which turns out to present significant deviations with 
respect to Eq. ()A4p . Indeed, if one takes for instance 
5i = 10 and (3 = 1.2, which are of the order of magni- 
tude of the experimental parameters in Ref. [3], there is 
about a factor two between the two expressions. 



APPENDIX B: BARGMAN REPRESENTATION 
ON THE TORUS 

Following [9], we present here the Bargman's repre- 
sentation on the extended torus. Classically, the two- 
dimensional toroidal phase space is a periodically re- 
peated cell having sides in suitable coordinates 
The classical dynamics is invariant under transla- 
tions by the elementary cell. Quantum mechanically, the 
states of the Hilbert space are required to be periodic 
functions (up to a phase) under translations defined on 
an extended torus, of size {MQ^P) (M is an arbitrary 
strictly positive integer). 



Tp\^) = e 



(Bl) 



where = (6>i,6>2) are two Bloch phases ranging from 
to 27r, and 



Tq = eyip{-iQp/h) 
fp = exp{iPq/h) . 



(B2) 



To simultaneously satisfy both equations ()Bip . Tq and 
Tp must commmute. This imposes that the area of the 
extended phase space, measured in units of Planck's con- 
stant /i, must be an integer 



MQP 

2iih 



N . 



(B3) 



The Hilbert space is thus an A/'-dimensional space 
parametrized by 6>, denoted Hn{0). For a fixed area 
MQP^ the semiclassical limit h is equivalent to 
N ^ oo. 

For each 1~Ln{0) one can define normalizable basis 
states \qn^O) and \pmi^) in the q and p representations. 
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respectively 



I, or Bargman representation, defined as 



\Pm,0) 

where 



V- 

Eoo 



-\qr, + vMQ) ,n = 0,...,Af- 1 
'■\p^ + vP) ,m = 0,...,Ar- 1 

(B4) 



N-1 



(B5) 



In Eq. (|B4p . \qn + vMQ) and \pm + I'P) are the usual 
position and momentum eigenstates, respectively, that 
have to be distinguished with respect to their periodized 
counterpart, \qn^O) and |pm,^)- The latter states sat- 
isfy the boundary conditions Eqs. ()Bip , and (n, 0\m^0) = 
ex.p{iqnPm/^)/ As Eqs. (|B5p show, the arbitrariness 
in the boundary conditions may be viewed as an arbi- 
trariness under shifs of order 1/N in the position of the 
discrete basis states \n^0) and \m^0) with respect to the 
intervals ^ : [0, MQ[ and p : [0, P[, respectively. 

Diagonalization of the Harper Hamiltonian in, for in- 
stance, the Iq^O) basis gives the eigenstates |'0i(6^)) = 
^n=o ^hri{0)\qn,0), z = 0, . . . , TV- 1, characterized by TV 

complex numbers that satisfy XI^Jq^ |'^i,n(^)P = 1- An 
alternative representation is in terms of coherent states 



(B6) 



n=0 



where 



{z\qn,0) = {TThy'/^exp (-1 [^{z^ + qD - V^zqr.]) 
e3{-'-i-i^iV2z-qn)[^) . 

(B7) 

The complex variable z = {q — ip)/ \pi denotes the central 
position of the coherent state, and 



+ CO 



6>3('u|r) = ^ exp(z7rrz/^ + 2i'i;z/) 



(B8) 



v= — oo 



is the Jacobi theta function [25]. 

Finally the Husimi representation Wr^^{z) of an eigen- 
state corresponds to the normalized square modulus of 
the coherent state representation 



(B9) 
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